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The impact of particle production during inflation on the primordial curvature perturbation spec- 
trum is investigated both analytically and numerically. We obtain an oscillatory behavior on small 
scales, while on large scales the spectrum is unaffected. The amplitude of the oscillations is pro- 
portional to the number of coupled fields, their mass, and the square of the coupling constant. 
The oscillations are due a discontinuity in the second time derivative of the inflaton, arising from a 
temporary violation of the slow-roll conditions. A similar effect on the power spectrum should be 
produced also in other inflationary models where the slow-roll conditions are temporarily violated. 

I. INTRODUCTION 

Inflation is considered one of the most promising candidates to explain the statistical features of the observable 
Universe revealed by the detection of the cosmic microwave background (CMB) @, [| and large scale structure 
surveys such as the Sloan Digital Sky Survey (SDSS) 

In the simplest models, inflation is an exponential expansion period of the Universe, driven by scalar field, 
called inflaton, slowly rolling down its potential. One of its main predictions is a primordial curvature pertur- 
bation spectrum of the form Pn(k) — k 71 "" 1 , where n s is the so-called scalar spectral index. 

Various extensions of the simplest models have been proposed, and in this paper we will consider the coupling 
of the inflaton to another massive scalar field, which has been previously investigated by different groups, leading 
to apparently conflicting results 0, Q. The first group was in fact obtaining just a local peak in the power 
spectrum for scales which leave the horizon around the time of particle production, while the second obtained 
a small scale oscillatory behavior leading to a step between large and small scales. We study the problem 
both analytically and numerically, obtaining an intermediate result which confirms the oscillation of the power 
spectrum on small scales, but without any step. We also derive an analytical approximation for the curvature 
power spectrum, which clearly shows the dependency of the amplitude and period of the oscillations on the 
mass and number of the coupled fields and the coupling constant. 

We mention that the presence of a temporary non-slow-roll stage during inflation and its effect on the scalar 
and tensor perturbation spectrum as well as the resulting CMB anisotropy have been investigated in @, 0| • Our 
analysis may be regarded as a Especial case where the perturbation spectrum can be studied analytically, and 
confirms other general studies [altJGil of the effects of singularities in the inflaton potential. See also a recent 
paper by Joy et al. [ll| for comparizon with the WMAP data. 

The paper is organized as follows. In Section HH we briefly describe the main features and motivations of the 
model we study. In Section [Ml we describe our analysis and present the numerical results for the spectrum of 
the curvature perturbation. In Section IIV1 wc adopt an analytical approximation and derive the spectrum in 
both large and small k limits. In Section [Vj we summarize the results obtained and provide some ideas about 
possible extensions. 



II. MODEL 



We will consider a theory with the potential, 

V(4>, p) = Vo + + ±N(m v - g<f>) 2 cp 2 , (I) 

where we assume Vq dominates over the other terms during the period of interest, so that the evolution of H 
can be neglected in the calculation. 

The equation for the inflaton can be written in the form Q : 

$ + mi + ml<t>-gN(m v -g<f>)(tf?) = 0; (2) 

(<p 2 )*9(t-t ) — ^-n W" 3 , no =g^ l -ff^-, (3) 
m v — g<p \ao J {Ziry 

where C is a constant of order unity, to is the time at which g<j> — m v , when the effective mass of the tp field 
becomes zero and most of the particles are produced. Here and in the following, the suffix denotes a quantity 



2 



evaluated at t = to. The two point function can be approximated with the expression above after renormalizing 
by subtraction of its asymptotic past value, when no particle were produced 0]. 

It can be shown that the time scale At c over which particle are produced is much smaller than the Hubble 
time, HA C ~ 10~ 3 « 1 5]. This justifies our approximation of (p 2 ) by a step function as given by Eq. for 
scales k/a < At^ 1 ~ 10 3 H. 

For later convenience, we slightly rewrite the above field equation as 



<t> + 3H<j> + ml4> = S(t) ; 



S = 9(t- t )gNCn — 
a 



We solve the equation for the curvature perturbation on co-moving hypersurfaces: 



z n ad> 

n'' + 2-n c + k 2 n c = 0; z = — 

z H 



(4) 
(5) 

(6) 



where a prime denotes the conformal time derivative, ' — d/drj — ad/dt. The inflaton perturbation on flat 
hypersurfaces 5d>f, which is to be quantized on subhorizon scales, is related to 1Z C as 



(7) 



We assume that there is no isocurvature perturbation after inflation, So, the power spectrum of the curvature 
perturbation is given by 



2wP% 2 (k) = V2tf\1l k (t f )\, 
where TZk is a properly normalized mode function, and tf is the time at which inflation ends. 



(8) 



III. CALCULATION 



First we solve the inflaton background. We assume the inflaton is slow-rolling before to- Introduce two 
independent homogeneous solutions of |4]), 



I7±(t)=exp[A±frt]; A ± =--(l Tv /l-V/9 

where /i 2 = rra^/iJ 2 , and the Wronskian, 

W{t) = U+U-- U- U+ = (A+ - X-)HU+U- = (A+ - \_)He~ 3Ht . 
the solution of our interest is expressed as 

/oo 
dt'G{t-t')S{t'), 
-oo 

where the Green function G is given by 

G = ^-*) wm ' 



<f>(t) = d, U+(t-t ) 

6(t-t )M 3 



(9) 



(10) 



(11) 



(12) 



Note that U+(i) describes the slow-roll solution, while U-(t) the rapidly decaying solution; A + w — /^ 2 /3 and 
A_ re -3 + m 2 /3 for fi 2 < 1. 

With the source term given by Eq. ([S]), we obtain the solution explicitly as 



A+A_(A+ - A_)# 2 



-\+U+{t ~ t Q ) + \-U-(t - t ) + (A+ - \-)e- 3H(t -^ 



(13) 



3 



where we have introduced the mass scale M by 



3/2 



M^gNCn^g^NC^- 

{2-k) 6 



(14) 



Taking the time derivative of this equation gives 



4> = \+H(f)QU + (t -t ) 
6{t - i )M 3 



A + A_(A + - X_)H 



-\\U+(t - to) + \ 2 -U-(t - t ) + (Ai - X 2 _)e- 3H (*"*«>) 



(15) 



We note that this result implies the presence of a step in </>. At late times, H(t — to) ^> 1, the above equation 
reduces to 



A+^o 1 



M 3 



A-(A+ - \-)H 2 <fto 



U+(t-t ). 



(16) 



Thus there is a step of the relative magnitude A0/0 ~ Af 3 /(9H 2 <fro) compared to the case of no particle 
production, and it will be reflected in the overall shape of the spectrum in general. However, as we shall see 
immediately below, for the values of the parameters we choose, the step turns out to be negligible. 
For ease of comparison, following [1,0, we will make the following choices for the model parameters: 



= 10 6 m p i , m ¥ 



2m 



pi ■ 



2™2 



.9 = 1, Vo = 5m m 



(17) 



where m p i = G x l 2 is the Planck mass. For this choice of the parameters, we find the step relative to the case 
of no particle production is small, A<fi/<j> ~ 10 _5 iV unless N is extremely large. The behaviors of <\> and <j> are 
plotted in Fig. 1 and 2, respectively. 
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FIG. 1: (j>(t) j (m p m^) is plotted for —0.1 < m^(i — to) < 0.5. The solid line corresponds to N = 1, the small dashed line 
to N = 8, and the long dashed line to N = 16. 



Setting u = aS(f>f — —z!Z c , where z = a<f>/H, we have 



u = 0. 



(18) 



Since we assumed that the potential is dominated by Vq, the time variation of H in z can be neglected and the 
scale factor a may be approximated by that of a pure de Sitter universe, a = (—Hrj)^ 1 . On the other hand, 
the time variation of <fi cannot be neglected, particularly at and after the transition. For </>, we use the solution 
given by Eq. (p~5|) with the identification Ht = \n(—Hr]). 
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FIG. 2: (j)(t) / (m p m^,) is plotted for —0.1 < mj>(t — to) < 0.5. The parameters are the same as Fig. 1. 



At rj = 770, 4> is discontinuous. Hence z" , which contains the third derivative <f>, contains a delta function. 
This implies v! is discontinuous at r\ = 7/0. To evaluate this discontinuity, we calculate the contribution of the 
delta function in z", 



— dr) = 

z 



a-o ,, n a Q M 3 a 
-r- = N gn C -y- = — , — 



(19) 



where </>q- and </>o+ are the values of 1 
at 77 = r]o for u is given by 



right before and after t = to, respectively. Thus the matching condition 



i o+ = u 'o- ivo) + D u - , u 0+ = U - 



(20) 



Turning back to the original variable 1Z C , it is noted that this matching condition implies that both 1Z' C and 1Z C 
are continuous at 77 = 770 - This is of course consistent with the evolution equation Q for 1Z C , in which there is 
no delta function. 

To calculate the power spectrum, we split it into two parts, corresponding to the modes greater or smaller 
than kg, where fco = (aH)o. We assume the standard Bunch-Davies vacuum for 8<pf at 77 — > —00, and solve for 
the positive frequency mode functions. Thus at sufficiently early times, 77 — s- — cxd, the mode function can be 
well approximated by 



p —ikri ■ 
~ (1- — ) 

V2k W 



(21) 



where it< denotes the mode function at r\ < rjo. For numerical analysis, we use this as the initial condition 
for each mode when it is inside the horizon and when 77 < t/q. Then we return to the original variable 1Z C and 
numerically integrate Eq. ^j. As we noted in the above, there is no delta function in this equation, but only a 
discontinuity in z' / z. Hence it can be numerically integrated across the time r/ — i]q without any problem. 

For modes k < ko, for which the sudden change in the effective potential of the inflaton happens after horizon 
crossing, Eq. (|21[) continues to be a good approximation until a mode crosses the horizon. Hence using it as the 
initial condition at horizon crossing, we solve the differential equation numerically. 

For modes k > fco, the particle production takes place before the modes leave the horizon. So, we set the 
initial condition at a sufficiently early time 77 — r\i < rjo common to all the modes, when the mode functions are 
well approximated by Eq. (|2ip . Then we integrate Eq. ^ numerically. 



The numerical results for the power spectrum are given in Figs. 3, 4 and 5. For the modes k < kg, we do 
not find any appreciable evolution of modes on super-horizon scales except for modes close to k — fco- For the 
modes k > k , the oscillatory feature of the spectrum is in agreement with [j| and other studies such as fl2l ]. 

Our result that there appears no step in the spectrum is in agreement with Adams et al. [HI , while it differs 
from Elgaroy et al. Q who claim that a step in the power spectrum is produced due to the non-conservation 
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^ (fc) is plotted for 3 x 10 < k/(aoHo) < 50 in the case of N = 8. The solid line is the numerical result, 



FIG. 3: R, 

the dashed line is the analytical approximation, and the long dashed line is the spectrum in the absence of particle 
production. 
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FIG. 4: The same as Fig. 3, but in the case of N = 16. 



of entropy perturbation. We disagree with this interpretation, since a source term should be present for the 
superhorizon evolution of 1Z C due to entropy perturbations, which is absent in the present case. 

In order to confirm our numerical results, we consider an analytical approximation to the problem in the next 
section. We find a good agreement between the numerical results and analytical approximations. 



IV. ANALYTICAL APPROXIMATION OF THE SPECTRUM 



In this section, we consider an analytical approximation for the power spectrum. As in the case of numerical 
analysis, we split the modes into the two, k < ko and k > ko, and discuss them separately. 

For the modes k < ko, since the transition occurs on superhorizon scales, the only possible time variation 
of 1Z C is due to the discontinuity in z'/z in Eq. ([B]). However, since z'/z — az/z, it grows exponentially large 
as time goes on and the discontinuity becomes totally irrelevant, unless the magnitude of the discontinuity is 
exponentially large (which is apparently not the case). Therefore, for k «C fcn , there can be no time evolution of 
1Z C on superhorizon scales. In fact, by using the technique developed in [Tj|], one can explicitly show that the 
spectrum can be modified only near k = ko and quickly approaches the standard result at k <C ko. Thus the 
curvature perturbation spectrum is given by the standard formula, 

P U\k) = -^— (k«k ), (22) 
2n\(j)(t k )\ 

where t k is the horizon crossing time, k — a(tk)H. 

As for the modes k > ko, the analysis is a bit more complicated. First we note that, a we stated before, the 
mode functions u are well approximated by Eq. (|21[) when r\ < rjo- After the transition, there is a short period 



6 



0.0006 
.00055 

.0005 
. 



A 








001 0.01 


. 1 


1 


1 10 



FIG. 5: P^ 2 {k) in the cases of both N = 8 and N = 16 are plotted for 3 x 10~ 2 < k/(a H ) < 50. Clearly the amplitude 
of oscillations is larger for larger N, approximately in proportion to N. The long dashed line is the spectrum in the 
absence of coupling. 



when the slow-rolling solution for the background inflaton does not hold. Nevertheless, for sufficiently large fe, 
for which the mode is still deep inside the horizon at the time of transition, this small violation of the slow-roll 
conditions is totally negligible because k 2 3> z"/z. Hence the approximate solution v defined in Eq. (f2"Tj) is 
still a valid solution even at ij > 770, including the non-slow-roll period right after the transition. However, v is 
no longer the positive frequency mode function there. Instead the desired mode function at 77 > 770 should be 
expressed as a linear combination of v and v*. Hence we may set 



= a k v + p k v* , 



(23) 



where u > denotes the mode function at r\ > r/Q. It is useful to note that the coefficients ak and f3 k can be 
expressed in terms of u>, v and v* as 



ak = — i(u*'it> 
This gives the formula for the spectrum, 



v w>) . 



H 2 



2n\<t>{t k )\ 



j3 k = i(u'w> - vuL) 



\otk - Pk\ 



It should be noted that \otk\ 2 

Now, the matching condition (f^TJ|) implies 



\Pk\ = !• This relation can be used as a consistency check. 



u >(Vo+) = v , u' > (r] 0+ ) = v' + D v . 



(24) 



(25) 



(26) 



Applying Eq. (|2T)) to an epoch right after the transition, rj = r/ 0+ , with u> and it^ given by these equations, 
we find 



Oik 



1 



iD v VQ = l + i 
D 



D 



1 



1 



-iDqVq 



2k 



2k V 1 (k Vo ) 2 
2 



1 



1 

kr]o) 



-2ikr)Q 



Inserting them into Eq. 

H 2 



we obtain the spectrum at k 3> ko as 



Pn(k) 



1/2 



2-K\4>{t k )\ 



£0 
k 



sin 2krjo + O - — 



(-)' 


n 2 r 




2k 2 _ 



1 + cos 2h] + O - — 



1 

h] 



1/2 



(27) 



(fc> fc ).(28) 



As seen from Figs. 3 and 4, the above analytical approximation agrees well with the numerical results at large 
k. However, the analytical approximation loses its accuracy for modes which leave the horizon around the time 
of particle production rjo- This is because the effect of non-slow-rolling is non-negligible for these modes and 
the approximate solution (|2ip is no longer valid. 

The spectrum at k > ko behaves like a dumped harmonic oscillator, where the amplitude of the oscillations 
is proportional to Do, provided Do/k < 1. As one can guess from the formula (|23|) . together with Eq. I|27p. 
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these oscillations are due to an interference between the positive and negative frequency mode functions, or 
contamination of negative frequency modes with positive frequency modes due to the transition. One might 
doubt if this is due to the approximation of replacing an otherwise smooth function by a step function in Eq. ([3]) . 
However, as we mentioned there, our approximation is valid for k < a^/ At = k${H At)^ 1 ~ 10 3 fco, while the 
oscillations are present for all fc > fco. This confirms that the oscillations are not an artifact of the approximation 
but real. 

Using the expression for Do, Eq. (fill)) , the amplitude is evaluated as 

£o _ OqAP_ _ a H Ng 2 m mj," 1 

k ~ k O _ k ^3(2^)3 HV 2 ' 1 ' 

This analytical estimate not only confirms the numerical result that the oscillations amplitude is proportional 
to the number of coupled fields N, but also shows that it is proportional to the square root of the mass m v 
and is proportional to the square of the coupling constant g. It also shows that the period of the oscillations is 
given by koir. 



V. CONCLUSION 



We have studied the impact of the coupling of the inflaton to a scalar field on the primordial curvature per- 
turbation spectrum. We found the presence of an oscillatory behavior on small scales, for the modes which leave 
the horizon after the time of particle production. We have also presented a very good analytical approximation 
for the evolution of small scale modes. This method can be applied to the general case of a sudden change in 
the inflaton potential, leading to a temporary violation of the slow roll conditions. 

The amplitude of the oscillations is proportional to the number of coupled fields N and the square root of 
their mass m^ 2 , and to the square of the coupling constant g 2 . The period of the oscillations is k^ir, where fco 
is the wavelength that crosses the horizon right at the time of the particle production. 

On large scales, k < fco, the power spectrum is virtually unaffected by the particle production, contrary 
to what was claimed in a previous numerical investigation [5J. This is because the violation of the slow-roll 
conditions can affect only on those modes close to k = fco on superhorizon scales, and the curvature perturbation 
is conserved on sufficiently large superhorizon scales, fc <C fco, no matter what occurs there. 

Our results are quite general in the sense that in models where slow-roll conditions are temporarily violated, 
the spectrum will have oscillations on scales smaller than the mode which leaves the horizon at time of transition, 
while it will remain unchanged on the larger scales. The presence of such features in the observed CMB spectrum 
could help to determine the magnitude and the lapse of periods during which the slow-roll conditions are violated, 
although it may be difficult in practice to distinguish such a feature of a primordial origin from a similar feature 
due to intermediate astrophysical processes happening before and/or after recombination. 
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